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ABSTRACT 



Within the simplest model of metals, namely a gas of electrons with Coulomb interactions, 
in the presence of a uniform background of positive charge to enforce electric neutrality of 
the system, we have derived a mechanism, by which the Coulomb interaction between the 
electrons generates a new kind of magnetism. The ground state of the metal is represented by 
a magnetically ordered state described by a non-local magnetic field. This non-local magnetic 
field does not produce spin polarisation of electrons, but induces a special long range correlation 
between electrons of opposite spin. This mechanism results in a theoretical value for the binding 
energy per electron, which is more than twice the corresponding value for the unmagnetised 
state of the metal. The new magnetic order proposed and analysed theoretically here, can in 
principle be experimentally tested. 
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1 Introduction 



In considering many body systems, especially condensed matter, using Quantum Field Theory, 
it is frequently advantageous to use functional integral methods rather than canonical operator 
methods. A particular benefit of functional integral methods is the possibility to eliminate 
certain fermionic functional integrals by means of a Grasmannian change of integration vari- 
ables, and to obtain in stead bosonic integrals involving new bosonic field variables. These 
bosonic fields can then be associated with existing macroscopic structures or phases appearing 
in interacting fermionic systems. To explain this in principle, consider a system of interact- 
ing fermions described by a given action Sf[tp* where ip*(r) and ip(r) are anticommuting 
(Grassmann-odd) variables, and r is a parameter in the interval {0,(3), with j3 the inverse 
temperature. 

All the termodynamical properties of the system in its equilibrium state are obtained from 
the partition function Z, which is expressed as a functional integral as follows, 

Z = J D(r^)ex P {-S f [r^}}, (1) 

where D{ip*,%p) is a standard measure over Grassmann variables [lj, @, ||. Making a change 
of integration variables in ([!]) known as the Hubbard- Stratonovich transformation, or more 
generally a change of variables which introduces Grassmann-even, i.e. bosonic variables B(t), 
one obtains a new functional integral representation 0, of the partition function (J!]) 

Z = J D(B,r,ip)exp{S bf [B,rM}- ( 2 ) 

The second representation (fj) is exact but not nesessarily useful or interesting, unless one 
manages to choose the bosonic field variables in some particularly convenient way. This depends 
naturally on the functional form of the original fermionic action Sf[ip*,ip}- There may exist 
several different actions Sbf [B, ip* , ip], corresponding to the same Sf[ip*,i/j). Under certain 
circumstances one may carry out the integration over the fermionic variables in (^) explicitly 
and exactly, obtaining an expression involving a bosonic field only, 

Z = jD(B)ex V {-S eff [B]}, (3) 

where S e ff[B] is an effective action functional involving only bosonic variables B. The functional 
integral ([|) can in general not be carried out explicitly, but may be very convenient as a a 
starting point for approximate calculations. 

We will presently use the machinery presented in general terms so far, to analyse the simplest 
model of metals, namely an electron gas with Coulomb interaction between electrons, in the 
presence of a uniform background of positive charge, which is used in order to enforce the electric 
neutrality of the system. However it is convenient to discuss the approximate evaluation of the 
functional integral (^) in general terms before considering detailed calculations. 

The approximate evaluation of ([3D goes in general as follows. One expands the action 
S e ff[B] as a power series in the fluctuations B = B — B around its stationary value S e ff[B]. 
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The mean field B is a solution to the variational Euler-Lagrange equations, 



5S eff [B] = 0. 



(4) 



Writing 



S eff [B] = S eff [B] + S fl [B] 



(5) 



one can then write, 



Z = Z < exp {-S f i[B}} > , Z = exp {-S eff [B\} 



(6) 



where the symbol < ... >o stands for average value with respect to S e ff[B\. In this method, 
the saddle point contribution is treated exactly, while the effects of the fluctuations around B, 
are treated perturbatively, to the desired accuracy. 

If the replacement of a fermionic integral with a bosonic one is succesfully carried out along 
the lines indicated above, then the bosonic fields B ought to describe macroscopic structures, 
or phases, of the interacting fermionic system, such as ferromagnetism, superconductivity and 
so on, for interacting electron systems. The effects of fluctuations 



ought to give rise to small perturbative corrections only, to the underlying structure of the 
system described by the field B. If one can ignore effects of fluctuations, then the partition 
function Z of the interacting fermion system is approximated by Z Q , and one says that the 
system is described in the mean field approximation represented by the field B. 

It was already mentioned above, that in general there may exist several realisations of 
fields B by means of which one transforms the functional integral ([!]) into the equivalent form 
(|2|). Hence, in the same fermionic system, one may in principle obtain several macroscopic 

structures, or phases, described by a set of mean filds B , where the index (i) enumerates 
the different macroscopic structures that can appear in the considered fermionic system, under 
given external conditions. With each mean field B one then associates a corresponding 
perturbation expansion, given by the relation @, for B = B^\ If there are several mean fields 
B , i — 1,2, ... one may ask which of these is actually realised in the system at hand. There 
is a nearly obvious answer to this; the structure B associated with the lowest energy under 
given external conditions is realised. It may of course happen that the external conditions do 
not lead to a unique structure B , in which case we have a coexistence of macroscopic phases. 

2 The simplest model of metals 

Let us now return to the physical problem of interest, namely the problem of the ground state 
in the simplest model of metals which we analyse using the general method outlined above. 
The problem in question is ancient ||, @, and concerns more precisely the question of whether 
the genuine ground state of the metal is an unmagnetised state or a ferromagnetic state. 



< exp{-S fl [B}} > 



(7) 
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The model in question consists of a large number N of electrons (an electron gas) in a finite 
but large volume V, interacting through the ordinary Coulomb interaction, in the presence of 
a uniform background of positive charge determined so that the system is electrically neutral. 
The system is described by the following Hamiltonian, 

H = H + H int — C, (8) 

where 



(9) 



H Q = [ ci 3 x£ 
describes the free electron gas, and 

Hint = £ E / / f^~4l {^(x')^(x)^(x)^(x')} (10) 
2 , Jv Jv x — x' 1 > 

describes the Coulomb interaction between the electrons. Here the quantities ip^ and ip a are the 
creation and annihilation operators, respectively, of electrons with spin component a. Finally, 
the constant C in the expression (^) is the following, 

C= e **[[^%, (11) 
2 Jv Jv | x — x' | 

where er] is the density of uniformly distributed positive charge related to the average density 
of the electron charge as follows, 

e V = eJ2<€(x)Mx) > ■ (12) 

a 

Incorporating the constant C as defined in ( |TTD above in the total Hamiltonian (§), makes the 
system electrically neutral, as it should be. 

The ground state energy £"0 of the unmagnetised state in the model above has been derived 



in perturbation theory ||, 0], ||, |10[ , with H as the unperturbed Hamiltonian , and expressed 
as a sum of three terms, 

corr. 

(13) 

The expression for the leading term Ekin. in (0) is as follows, 

e 2 N 3 /9tt\§ 1 _ e 2 N 2.21 



where is the total number of electrons in the system, a® is the Bohr radius, r s = r /a and 
r is a parameter related to the volume per electron, i.e. |vrrQ = V/N. 



The term E exch . in ( ]T3] ) is the exchange energy, i.e. the expectation value of H in t. in the 
ground state of a free electron gas, and has the following form, 

e 2 N 3 /9tt\ i 1 e 2 N 0.916 . . 

Eexch. = -15— ir\~r) ~ ~ ~^ • 15 

2a 27r V 4 / r s 2a r s 
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The last term in the expression ([13]), is the so-called correlation energy, a term introduced by 
Wigner 0. This term was calculated by Gell-Mann and Brueckner to the accuracy e 4 , by 
summing an infinite number of so-called ring diagrams, with the result, 



e 2 N 
2a n 



2 



— (1- log 2) log r s - 0.096 



7T 2 



e N 
2a 



[0.0622 log r s - 0.096]. (16) 



It was pointed out by F. Bloch already in 1929 that the ferromagnetic state of the electron 
plasma in metals could have a lower energy than that of the unmagnetised state. By considering 
a polarised gas of electrons with N + electrons in the spin-up state, and N_ electrons in the 
spin-down state and defining the polarisation £ = (iV + — N_)/N, where N = N + + iV_, one 
finds the following expression for the ground state energy to the first order in the interaction 
potential, 

e 2 N 3 /9n\ 



- 



2a 10 V 4 J 
e 2 N 3 /97T\^ 



;i+£)3 + (i-o 



2a An V 4 



;i + o^ + (i-o 
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One readily infers from the expression (|I1]) that the ferromagnetic state, with £ = 1, has a 
lower energy than the unmagnetised state with £ = 0, provided that 

2tt /9tt\ 3 r i. . . 

r s > — [-) [23 + 1] « 5.45. (18) 

This problem is also discussed in the textbook by Fetter and Walecka |]]. 

In the relation ( |i~7D the polarisation £ is a free parameter subject to variation in the interval 
[0,1]. The expression (|TT| ) represents, in fact, only the energy of the electrons in a uniform 
magnetic field, which gives rise to the polarisation of these electrons. If the electrons become 
spontaneously polarised, then there should appear a spontaneously self-organised magnetic field 
H in the system, which makes a positive contribution Ef^d to the total energy of the system. 
Then the total energy of the ferromagnetic state Ef erro would be the following, 

Ef err o E ma g -\- Eji e id (19) 

There is the possibility, that the energy Ef erro may be lower than the energy of the unmagnetised 
state given to lowest order by the sum of (0) and fllSP . Thus the ferromagnetic state can in 
principle be a genuine ground state of the system for certain external conditions. 



In m and |10| doubts were expressed concerning the possibility of having a ferromagnetic 
state as a genuine ground state of the Hamiltonian (|S[) in view of the fact, among other things, 
that the contribution Ef ie i d was neglected in the considerations based on ([H 



Independently of arguments for or against the possibility of a ferromagnetic ground state 
for the system described by the Hamiltonian (||), the Bloch considerations based on (|l9|) pose a 



4 



challenge in electron plasma theory. The challenging task is to derive a mechanism by means of 
which the Coulomb interaction between electrons in the model described by (|8]), can generate 
magnetism. To the best of our knowledge, this has been an unsolved problem so far. Below we 
give a solution to this problem, i.e. we show, using the functional methods discussed in general 
terms in the Introduction, that there is a kind of magnetic field that arises due to the Coulomb 
interaction, and that the corresponding magnetic state under certain conditions is a genuine 
ground state of the system under consideration. 



3 Magnetic ordering in the electron plasma 



The grand canonical partition function Z corresponding to the Hamiltonian (H) is given by the 
functional integral (0) with the following action 



S/ [</-*,*/>] = 

,2 r p 



(20) 



d h 2 

C( x > r)—i/; a {x, t) - 7^C( X > t)V 2 ^ q (x, r) - ^(x, r)^ a (x, r) 
d 3 xcf 3 x' 



+ \J / 1 x "| [^(x / ,rX(x,T)^ a (x,T)^(x / ,T)] - 0C, 

® a,a' I I 

where the electron field variables ^>*(x, r) and t/> q (x, r) are Grassmann variables enumerated 
by the spin index a, the space coordinate x and the "time" r, and fi is the chemical potential. 



The electron field variables satisfy the antiperiodic boundary conditions, 
4>* a (x, t + (3) = -C( x > T ) > ^(x, r + (3) = --0a(x, r). 



(21) 



3.1 The unmagnetised electron plasma 

We first review one known , [0] way of introducing bosonic variables in the fermionic system 
described in terms of the functional integral ([!]). This consists of introducing an electric scalar 
potential 0(x, r) , in the present case when the action is given by the expression (BO) above. 

As mentioned in the Introduction, there may in general be more than one way of transform- 
ing the functional integral ([]]) into a functional integral (§) involving also bosonic (Grassmann 
even) variables. We consider an alternative to the procedure explained below in the next 
subsection. 

The partition function Z involving the electric scalar potential is the following, 

Z = ±- [ Dtfrf^expi-SvfaPM, (22) 
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where 



S bf [(p,ip*,ip} 



±- f dr I rf 3 x (V0(x,r)) 2 
Sir Jo J 

+ J o drj d 3 x ^^(x,r)- ^ a (x,r) 

+ / dr I d 3 xV 
Jo J 

+ie J dr J d 3 x 0(x, r) 



(23) 



2m 



In Eq. ([23|) above, the quantity erj is a uniform positive charge density and Zq^ is defined by 
the following functional integral, 



Z H = J L»(0)exp|-i-^dT | rf 3 x [V0(x,r)] 



(24) 



The quantiy 0(x, r) introduced above, is a bosonic, i.e. Grassmann even, integration variable, 
which satisfies the following periodic boundary conditions, 



(x, r + P) = 0(x, r) 



(25) 



In the action fl2"3|) one can recognize a contribution from the energy density of an electric field 
E(x, r) = — V0(x, r), and from a system of electrons interacting through a scalar potential 
0(x, r). Performing the integration over the variable 0(x, r) in the functional integral ( p2|) one 
recovers indeed the functional integral ([I]) with the action S given by Eq. (|2T|). However, one 
may also first integrate over the variables ip* and ip in Eq. (E3). One then obtains, 



1 



-D [0] exp 







dr / rf 3 x 



1 



— [V0(x, t)] + ie?70(x, r) 



■57T 



det ^(0) I , 



(26) 



wherere ^"(0) is a functional matrix, the matrix elements of which are as follows, 

d h 2 



< x , r , a' | J r (</>) | x, r, a > 



dr 2m V ^- ie0(X ' r) 



5(r' - r) 5(x' - x) <W (27) 



The expression (|26|) defines a certain effective action 5" e //. Evaluating the functional integral 
(|26"D approximatively, to the order e 4 , using the general procedure given by Eqns. (U) and (§) 
above, one reproduces the expressions fll~4]), (|I5]) and ( |T6"D in the expression ( |T3| ) for the ground 
state energy of the unmagnetised electron plasma. 



3.2 The magnetised electron plasma 

We now come to the central issue in this paper, namely that there is another way of recognising 
bosonic variables in the system described by the action fl20|) , than the one described above, 
which utilised a scalar field </>(x, r). The novel transformation involves bilocal fields. 
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Consider auxiliary commuting (Grassmann even) real variables A a > a (x!, x, r), enumerated 
by spin indices a, a' = ±, by two space coordinates x' and x (x' ^ x) and by a "time- variable" 
t. The commuting variables A a i a (x.', x, r) are taken to be antisymmetric with respect to a 
simultaneous transposition of space- and spin indices, 

A*'a(x', x, r) = -A aa , (x, x', r) (28) 

Furthermore, these variables are supposed to satisfy periodic boundary conditions, 

A a / a (x / ,x,r + /3) = 4,» a (x',x,T) (29) 

The variables A a / a can be regarded as matrix elements of a certain 2x2 matrix, expressed 
with the aid of the 2x2 unit matrix a = 1 and the usual Pauli matrices a k , (k = 1, 2, 3), and 
a four-component quantity B = (B , B 1} B 2 , B 3 ) as follows, 

3 

A a , a (x', x, r) = *i>aBe(x', x, r) (30) 

£=0 



The symmetry conditions ( p28|) satisfied by the real quantities A a i a imply the following reality 
and symmetry conditions for the four- component quantity B, 

Im B k (x. f , x, r) = 0, B k (x.', x, r) = -B k (jc, x , r) ,k = 0,1, 3. (31) 

The second component B 2 is purely imaginary, and symmetric under the interchange of x and 

Re £ 2 (x', x, r) ee 0, 5 2 (x', x, r) = + £ 2 (x, x', r) (32) 
The real variables i a ' a (x',x, r) are used to define an action So& as follows, 

S ob [B] = -Lfdrj J d 3 xd 3 x' | x - x' | E 4k« ^ t) (33) 



x' 



= -^fdrff d 3 xd 3 x' | x - x' | E 4(x', x, r)5 fe (x', x, r), 

fc=0 

and the corresponding partition function, 

Z ob = J D(B)exp{-S ob [B}}. (34) 

We wish to evaluate the partition function (H) with the action Sf given by Eq. (p0|) . To this 
end we multiply Eq. ([I]) with the evident identity, 

1 = -L / D^exp^S^fl]}. (35) 

This results in the following expression for the partition function Z, 

Z=±- [ D(B)D(r,i>)e W {-S 0b [B] - SfWM (36) 



a' a 
3 



In the expression fl3"6|) we finally make the following change of integration variables, 



A Q / Q (x',x,r) -> A^ Q (x',x,r) = A,, a (x',x,r) + (37) 

#a'( x ', 7")V>a(x, r) - V* (x, r)Vv(x', r)]. 



ie 2 



The reason for the shift of integration variable described above, is simply that by evaluating 
the functional integral (^) in the manner indicated, one cancels the Coulomb interaction term 
in the expression (|20|) for Sf[ip*,ip] and is left with the following simple expression, 



Z=^- [D(B,r,ip)exp{-S bf [B,r,i>}} (3* 



with an action S b f[B,ip* ,ip] given by the following expression, 



S bf [B,r^] = S ob [B}+ |rf 3 x^*(x,r 



;2 



<9r 2m ^ 



^(x,r) (39) 



+*jf/ /rf 3 xd 3 xV*(x',r) ^t7*S fc (x',x,r)^ V>(x,r)-/3C 



In the expression fl3~9|) we have used ordinary two-component notation, 



for the Grassmann variables t/>*(x, r) and ^>(x, r ) ; respectively. 



V>= ( ^ + ) , = (40) 



The action (|39f) apparently describes a system of electrons interacting through a bilocal 
scalar potential field -Bo( x ', x, r) and a bilocal magnetic field B(x', x, r). The special symmetry 
properties of these field variables under interchange of space coordinates are given in Eqns. 
SID and (m). 



The formal transformations ( |3"S| ) and ( |3"TD leading from the original action ( |2"0"| ) to the ac 



tion (|39|), which describes a system of electrons interacting through a bilocal field -B(x, x',r), 
define a mechanism through which magnetic interactions can be generated from the Coulomb 
interaction. This mechanism is an alternative to the transformation involving only an electric 
scalar potential 0(x, r) which resulted in the action S b f given in Eq. (|23|) , by means of which 
one could reproduce the expressions (|14"D, (|T5|) and (plj) in the formula (|T3|) for the ground state 
energy E of the unmagnetised electron plasma. Clearly, the self-consistency of the mechanism 
described above will have to be ascertained. We return to this question presently, and continue 
now with a straightforward analysis based on the action ([39|), which is bilinear in the fields ip* 
and ip. This circumstance makes it possible to carry out the functional integrations over ip* 
and i/j in the expression (|38"D exactly and explicitly, with the result, 



e 0C 



■'Ob 



j D[B] exp {S ob [B}} det F{B), (41) 
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where F(B) is a functional matrix, the matrix elements of which can be written as follows 

d h 2 



< x', t', a' | J-{B) | x, r, a > 



dr 2m ^ ^ 



5a>a S(t' - r) <5(x' - x) (42) 
+i 5>i, a ^(x',x,T) 5(t'-t) 



1=0 



The Dirac 5-function 5{r' — r) which enters in the relation (f42f) above, is here defined in terms 
of its spectral representation, 



1 +°° 7f 

8(t'-t) = - exp^(r'-r), lq v = ~{2v + l) 

P u=-oo P 



(43) 



where the summation index v runs through all integer values, and the quantities u v are the 
so-called Matsubara frequencies. 



Employing the relation 

det F(B) = exp{Tr log F(B)} 
we obtain the final result for the partition function Z in Eq. (|41| 



Z=-L [D(B)exp{-S eff [B]} 



(44) 



(45) 



where the effective action S e ff[B] is as follows, 



S eff [B) = S ob [B] ~Y. dT [ rf3x < x, r, a | log F{B) \ x, r, a > -pC. 



(46) 



The equation (|46|) will be the starting point for an approximate evaluation of the partition 
function Z under consideration. This problem will be considered in the next subsection. 



3.3 The saddle point equations 



We now look for the extremals of the effective action ([4"6]), i.e. for solutions to the equations, 

5S eff [B] = 0. (47) 
It is reasonably simple to evaluate the variations involved in Eq. (|47|). Taking in addition into 



account the reality and symmetry conditions ([31]) and (|32D above, one obtains the following 
four equations from the conditions ([47]). For £ = 0, 1, 3 the equations in question are, 



— I x — x' I i?/(x, x', t) 



(4£ 



- ^ o'aa' { < x ') r > oc' | F~ l {B) | x, t, a > — < x, t, a' | jF _1 (i?) I x', r, a > \ 
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For I = 2 one gets 
1 



x — x' I .Bofx, x', t) 



(49) 



— - ^ a aa> {< x '> r,a' \ J 7 1 (B) | x, r, a > + < x, r, a' | T 1 {B) | x ; , r, a >} 



Using the relation ( ^2|) one finds readily that the quantity < x', r', a' \ T 1 (B) | x", t", a" > 
satisfies the following equation, 

J2 J rf 3 xO(x, a; x', a') < x', t', a' \ T~\B) | x", r", a" >= 5(r' - r> 3 (x - x")<W' (50) 



where the operator O is the following, 

d h 2 



(9(x, a; x', a' 



9r' 2m 



<W5 3 (x - x') + % °L B e(x, t')\ . 

e=o ) 



(51) 



The equations (fllf ) and (^) together with the equation (|50| ) constitute a complicated set 
of nonlinear integral equations for the saddle point configuration -B^x, x', r), \i = 0, 1, 2, 3. It 
remains to show that a nontrivial saddle point configuration exists. We shall not attempt here 
to discuss the most general solution possible of the set of equations fl48|) , (f49"D and ([5(]). Rather 
we confine our attention to the following special r-independent ansatz, which we denote by B, 

5,(x, x', r) = 0, i = 0, 1, 3; 5 2 (x, x', r) = -itffl x - x' |), (52) 

where H(\ x |) is a real- valued function to be determined. Inserting the ansatz (|52|) in the 



equation (|50|), one finds that this equation, apart from a factor —1, is nothing but the equation 
for the Green's function for an electron in a magnetic field of variable magnitude H in the 
direction of the 2-axis. Under these circumstances it is appropriate to introduce the following 
notation, 

G a i a " (x' - x", t' - t"; H) = - < x', r', a' | .F" 1 (B) | x", r", a" > . (53) 



We consider first the equation (f49|), which becomes an equation for the determination of 
the function H (| x |) in the ansatz (|5*2"D. In order to analyse this equation, it is convenient to 
consider the Fourier transform G aa i of the Green's function G aa > defined by the equation (|53"D, 

G aa / (k, v\ H) ee [ dr [ d 3 x G QQ ,/(x, r; if) exp(— zo^r + ik • x). 
jo iv 

Using the equation ( |50"D in the present case, one readily finds that 

f ,h 2 k 2 > 



E 



2m 



Li)]5 aa > - a aa ,H(k) \ Gv a »(k, z/; if) = 5 Qa », 



where 



H(k) = / d 3 xii(| x |) exp(zk • x). 



(54) 

(55) 
(56) 
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The Eq. (|55|) has the following solution, 

G aa , (k, v\ H) = - — — — -, (57) 

(zu-v - £+(k) + /^(mv - £-(k) + /i) 

where 

fi 2 k 2 fi 2 k 2 
S (k) = " £ ± (k) = — ± F(k). (58) 
2m 2m 



In passing we note the symmetry of the Green's function in Eq. ( |53| ) under the interchange 
of x' and x", which readily follows from the equation (|57|) above. From this symmetry follows 
that the ansatz (|52|) is consistent with the remaining equations (fP|). We then continue with 
our analysis of Eq. (E9[) . 



The equation (|49| ) now finally boils down to the following 

1 



e 2 



x | H(\ x |) = -E^^ a (x,0;if), (59) 



where 

1 * ~ 

G a / a (x,0;if) = — V lim V G Q ' a (k, z/; H) exp(-ik • x). (60) 
Inserting the expression ( (5T| ) in Eq. ( |60[ ) and using the well-known series, 

N 1 11 

lim V = (61) 

^ ^iv Z7r ( 2z/ + !) - ^ E exp(/5E) + 12' K ' 

one obtains from Eq. ([59]) the equation for the unknown function H in its final form, 

x I H(\ x |) = t £ { — ^L— ^ - — expHk ■ x). (62) 

The equation (|62"D is still a fairly complicated nonlinear integral equation for the determination 
of the function H; it should be remembered that the if -dependence on the right hand side of 
Eq. ( |62"D occurs through the energy spectra E± defined in Eq. (|58]) . In the next subsection we 
consider a limiting case of Eq. (0), for which the solution H can be constructed explicitly. 



3.4 Solution of the integral equation at T = 

In this sub-section we consider a limting case, namely the limit T — > 0, of the integral equation 
(p2]). This limiting case, which corresponds to (3 — > oo, is physically relevant, and furthermore 
such that one is able to obtain a fairly complete picture of the solution to Eq. (|62|) in the limit 
in question. 
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Let us first note that Eq. ( |5^) is equivalent to the following equation, 
I A | x | H { \ x |) exp( !k • x) = a* _ __L__ } . (63 ) 



In the limit (3 — > oo we have, 



lim r = S(£ F - E) (64) 

where G is the ordinary step function, and we use the notation, 

lim u = £ F (65) 

In the limit T — > the analysis of the equation (^) boils down to finding the roots of the 
following equations, 

£±(k)=£ F (66) 
where the quantities £±(k) have been defined in Eq. (|58|). 

It turns out that the relevant parameter in the analysis of the Eqns. (^6|) is the value of 
the quantity H(0), which, without essential loss of generality, can be taken to be non- negative, 
as is easily verified. One still has to consider separately the cases when £p > and Sp < 0, 
repectively. The necessity to consider also the case with £p < is a consequence of the fact 
that a part of the quasiparticle spectrum £_(k) is negative. We consider first the case 

£f>0 (67) 

There are then two cases to be considered, namely 

< H(0) < £p (68) 

and 

H(0) > £p (69) 



We first consider the case when Eq. (|68"D is in force. Then, assuming reasonable smoothness 
of the function H(k), both the equations ( |66D can have roots k = k < and k = k>, respectively, 
where 

h 2 

£ + (k<) = — k 2 < + H(k < )=£ F (70) 
2m 

and 

E4k > ) = £-k£-H(k > )=£ F (71) 
2m 

In the limit (3 — > oo one then gets the following equation from Eq. (|63|) , 

\ drr 2 H(r)sm(kr) = e 2 —{e(k > -k)-e(k < -k)} (72) 
jo Ait 
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From Eq. (|72| ) follows readily that 



H(r) 



2vr 2 r 4 



{[sin(fc>r) — fc>r cos(/c>r)] — [sin(/c<r) — k<r cos(A; < r)]} , 



(73) 



from which follows, 



H(k) 



e 

71 



2k 



k 2 

log 



k> + k 



k 2 -k 2 , 



k K + k 



(74) 



It then remains to determine the actual values of the roots /c< and k> from the equations ( |70D 
and (|71~D, respectively. From Eq. ([74]) one readily gets the following results, 



and 



iJ(*<) 



e 

7T 



e 

7T 



2fc. 



2A; 



k 2 

— log 



k 2 

— log 



/l> -|- /c< 



An> -|- /c< 



(75) 



(76) 



Inserting the values for H{k < ) and H{k > ) from Eqns. (I75p and flTB| ) into the Eqns. ( |70D and 
(|7T|), respectively, one finally obtains a pair of transcendental equations for the determination 
of the roots k < and /c>. The analysis of these equations is facilitated by the introduction of a 
momentum variable k + , related to the Fermi energy Sp as follows, 



h 2 . 
— / 

2m 



(77) 



It turns out that the existence of solutions to the transcendental equations ( [70]) and fl7l[) 
determining the roots k< and k>, depends on the values of the following dimensionless parameter 
aio, 

a = —V (78) 



where a denotes the Bohr radius, 



a 



7ia k A 
me 2 



(79) 



It is not difficult to verify that the transcendental equations ( |7D| ) and ( |7T| ) have roots k < and 
fc>, with k y > k< > 0, if and only if the parameter ao introduced in Eq. (|78|) satisfies the 
following inequality, 

\/6 

< a < « 0.204. (80) 

The inequality (|80|) yields a rather large lower bound on the Fermi energy £p. Namely, ex- 
pressing the scale k + in the expression (|77|) in terms of the parameter a defined in Eq. ([78]), 
and using the inequality ( |80D we obtain, 



S — — — — 

7t 2 «q 2a 7r 2 2a 



^1) 
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The lower bound ( pip is numerically far above the Fermi energies for all existing metals, 
completeness we nevertheless continue with the analysis of the case at hand. 



For 



The pair of transcendental equations (|70|) and (71) can only be solved numerically. For this 
purpose we introduce dimensionless parameters u and w by writing 



uk+ , An 



12) 



In Table 1 below, we give values of the parameters u and w introduced above in (|82"D , which have 
been obtained by solving Eqns. ( |70|) and ([71]) numerically, for a set of values of the parameter 



ao in the interval (80). 



Table 1: k< = uk + , k > = wk + 





a = 0.1 


a = 0.15 


a = 0.2 


u 
w 


0.98120 
1.01794 


0.88503 
1.08974 


0.61358 
1.19723 



In view of the condition (|3~1~D, the solution discussed above is presumably irrelevant for 
realistic physical systems. Hence we turn our attention to the case when the inequality (|69| ) is 
in force, which is actually a much simpler case than that already considered. 



We now assume that the inequality ( |5"9"D is in force. Then, still assuming reasonable smooth- 
ness of the function -ff(k), one finds that Eq. (|70|) has no non-negative root at all. Hence the 
equation (B3I) can now be written as follows, in the limit T — > oo, 



k 

drr 2 H(r) sin(fcr) = e 2 — 0(fc> — k) 



13) 



where fc> is determined by the equation (|7l|) alone. From Eq. (|83D one readily obtains the 
following expression, 



H(r) 



2vr 2 r 4 



{sin(/c>r) — /c>r cos(A;>r)} . 



which immediately leads to the following expression for the function H(k), 



H{k) 



g 2 r £2 _ £.2 

7 *> + T log 



k> + k 



We then finally obtain 



and 



H(0) = lim H(k) 



2e 2 



7T 



Hfa) = - fc>. 

7T 



(84) 

(85) 

(86) 
(87) 
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Using the result fl8~T|) in Eq. (|7l|) one obtains a simple quadratic equation for the determination 
of the root k > . It is still expedient to use the momentum scale k + defined by Eq. ( |77D as well 
as the dimensionless parameter a defined by Eq. (f78[). Writing 



k> = wk+ (88) 
one finds that Eq. (|71|) is equivalent to the following equation, 

w 2 - 2a w -1 = 0, (89) 
which gives the final expression for the root &;>, 



k> = K + V«o + 1) k+ = (1 + VI + «o ) (9°) 

v 7ra 

Let us recall that the analysis above is done under the assumption that both the inequalities 
( 1571 ) and (^) are in force. The latter inequality is now equivalent to the following condition on 
the dimensionless parameter a , 

«o > ^ « 0.204 (91) 

For future reference we note that the inequality ( |9"ID implies that 

2 fi 
— <k > <— (92) 
7rao vrao 

It remains to consider the case of negative Fermi energy Sp, 

Ep < (93) 

In this case, there is only one equation for the determination of the root k>, namely Eq. fl7lf). 
This equation can have a positive root fc> under the condition fl93"|) , if and only if 

H(0) > -E F (94) 

The analysis proceeds along the same lines as in the previous case with only minor modifications. 
We introduce a momentum scale k_, determined by the (negative) Fermi energy, 

h 2 

E F = -Tz-kt (95) 
2m 

as well as a dimensionless parameter 70 analogous to the parameter ao in the previous case, 

7o = ~\- (96) 

Writing 

k> = vk. (97) 
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we now obtain the following quadratic equation from Eq. ( |7T|) for the determination of the 
parameter v 

v 2 - 2 l0 v + 1=0 (98) 

This gives two roots k 



K = (7o - f^l) k- = (l- Vl-70" 2 )— (99) 

and 

K = (7o + VTo 7 ^) k- = {l + ^jl-% 2 ) — (100) 

7rao 

There is the obvious requirement that the roots k> be real; from this follows that we must 
necessarily have 

7 2 > 1 (101) 
The equations (0), ( """"D and ([8*6]) are now valid as such, but with fc> replaced with either the 
root given in Eq. (|99|) or in Eq. (|100|) . For negative Fermi energy £p < there are thus 
two distinct solutions to the equation ( |6"3"|) in the limit T — > 0, for each fixed value of the free 
parameter 70 satisfying condition (|101|) . 



The inequality (|9J) does not give any additional restrictions on the parameter 70 besides 
the restriction (101). For future reference we finally note that the roots fc> satisfy the following 



inequalities, as a consequence of the inequality ( |101| ), 

< k~ < — (102) 
na 

and 

< k+ < (103) 

7rao vrao 

Let us summarise the situation so far. We have shown the existence of solutions H(r) to the 
equation ( |62] ) in the limit T — > 0. If the Fermi energy £p is positive, then the solution iJ(r) takes 
the form given in Eq. (|73| ) if the inequality (168~|) is in force. This is the case if the dimensionless 
parameter ao defined by Eq. ( [78] ) satisfies the inequality (P0[). If the inequality ( "jlaf ) is in force, 
then the solution takes the form given in Eq. ([84]). This is the case if the parameter ao satisfies 
the inequality (OT). The latter solution can be obtained from the former in the formal limit 
k < 0. In both cases the solutions depend on a free parameter ao constrained either by the 
inequality (|8~0| ) or by the inequality Q9l|) . 



If the Fermi energy £p is negative, then the solutions H(r) take the form given in Eq. (p4|), 
with either /c~ given by Eq. ( p9|) replacing the parameter fc>, or /c+ given by Eq. ( |100| ) replacing 
the parameter k>. 

The demonstration given above, is by no means a replacement of a fullfledged proof of 
existence of solutions H(r) to the equation (|62]) in the general case, i.e. for nonzero finite T. 
However, the demonstration above is a proof of existence of solutions H(r) in the limit T — > 0, 
by explicit construction. The solutions obtained will be used below to calculate the ground 
state energy of the system under consideration. Before that, we consider in the next section 
the effects of fluctuations around the solutions discussed in this subsection. 
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4 Fluctuations 



We have found a stationary point B of the effective action (f46[) , namely the ansatz (|52D, in 
which the function H is the solution to the equation fl62|). We then have to consider the effect 
of fluctuations around this stationary point. We define a fluctuation field B as follows, 

fi(x, x, r) = B(x, x', r) - B(x, x', r) (104) 



Then one has to expand the effective action fl4"6|) in powers of the field B. We will make this 
expansion up to second order in the field B. After some calculations one obtains the dollowing 
result, 

S eff [B] = S eff [B] + S ob [B] + Sfi[B] - PC (105) 

where 

Sfi[B] = f dr fdr'f d 3 x J rfV J d 3 x"| dV' £ £ (106) 

The quantity in the expression above is given in terms of the Greens functions G defined 
in Eq. (EBI) as follows, 



R[B] = - l -G aa , (x - x', r - r'; B M (x', x", r')G^ (x" - x"', r' - r; H)a v a „, a B v {x!" , x, r) 

(107) 

The effective action S e ff[B] at the stationary point B has the following explicit form, 

S eff [B] = A, f d 3 xd 3 x' | x - x' | # 2 (x - x') (108) 



2e 2 

- ^log{[l + e- /3 ^«-^][l + e-^«^]} 



Let us now consider the physical meaning of the stationary point B in some detail. The 
effective action S e ff{B] originates from the following functional integral, 

exp{-S eff [B]} =exp{-S ob [B]} J Vty*, exp {-S of [i>*, 1>)} (109) 
where S f [ip*, ip] is the following action, 

S of [r,ip] = / o /3 rfr|rf 3 x|^(x,r)^(x,r)+^(x,r) 

+ J q drj d 3 x J rf 3 x / iJ(|x-x , |)^(x , ,r) ( T 2 ^(x,r) 

The action ( [L10|) describes a system of electrons interacting with a given non-local magnetic 
field H = (0, H(\x — x'|), 0). This magnetic field generates a correlated creation of an electron 
with spin up (down) at the point x', together with a simultaneus annihilation of an alectron of 
opposite spin down (up) at the point x. 



2m 



^(x,r) (110) 
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Introducing the following parametrisations for the Grassman variables ip and ?/>*, respec- 
tively, in Eq. (g0|), 



^(x, r) = -=Y: «k(r)e 4k ' x , ^(x, T ) = —J2 a^e^ {III) 
V K k vv k 

one can express the action ( |110| ) in terms of the variables a k (r) and ^(t) as follows, 



jf *"E {<(^)|: a k(T) + a*(r)[£ k - /i]a k (r)| (112) 



+ / dT^#(k)a k (r)(7 2 a k (T) 



r/3 

dr 

k 

where if(k) is the Fourier transform of the function if (|x|) defined as in Eq. (|56|) . 



The Grassmann variables a k a( r ) and a kQ ,(r) are associated with electron states specified by 
the wave number k and spin a = ±. In order to diagonalize the action ( |112j ) we introduce new 
fields c kce and c ka by the following unitary transformation, 

1 1 
c k+ = "^/|( a k+ - «Qk-) , c k _ = -^=(-m k+ + a k _) (113) 



"k+ " y^(° k + * fl k-) 5 c lc- — /^{^ a h+ + a k- 



Then the action ( |112j ) can be expressed as follows, 

r/3 



S of (c*,c) = J dr^U 



d 



c k + £ + (k)c* + c k+ + £_(k)c k _c k _ (114) 



where the quantities £±(k) are the energy spectra (|58|) of the quasiparticles described by the 
fields ( |113| ). As evident from the relations ( |113| ), the states of these quasiparticles are special 
superpositions of electron states with spin up and spin down. The probability to find an electron 
in a spin up or down state is equal to \ for each quasiparticle spin state. Finally, carrying out 
the appropriate functional integration in Eq. (|109|) using the form (|114|) of the action, one 
obtains 

exp {-S eff [B}} = exp {-S ob [B\} J] {[1 + e -^+«^][l + e -/^-(k)-M)]} (115) 

k 

which correponds exactly to the result given in Eq. (|108|). 



It behooves us still to check whether the self-organised magnetic field H = (0, i/(|x — x'|), 0), 
which is of order e 2 , survives under fluctuations which are also of order e 2 , or not. This involves 
an evaluation of the functional integral (|4"5|), with the effective action given by the second order 
approximation (|105|) . One obtains the following result, 

Z = exp{-S eff \B\}z fi (116) 
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where S e ff[B] is given by Eq. ( pL08| ) and 

[3e 2 rd 3 xd 3 x' 



log Z ft = ^ J P>G(x ~ x > °+; H )Y 



117) 



In Eq. ( 117|) above, the notation Tr stands for the appropriate spin index sums of the indicated 
quantities involving the Green's functions defined by Eq. (|53|). The term f3C in Eq. (|105|) was 
exactly cancelled by the direct Coulomb term which appeared in the average value < Sfi(B) >q. 



Using the results (|116|) and (|117|) on gets the grand canonical potential Q, 

Q = —ksT log Z 



1181 



from which one obtains the energy E mag of the electron system with the self-organized magnetic 
field H, 



E 



■may 



1 

2^ 



cf xdV I x - x' I H 2 fx 



(119) 



+ E {Mk)/(Mk)) + f_(k)/(f_(k))} 



- ^^!^ [TrG(x _ x , )0+;i/)] 2 



where the function f(S) denotes the Fermi function, 

1 



e m-p) + 1 



(120) 



In the expression ( |1 19| ) the first term corresponds to the energy of the self-organised magnetic 
field H,and the second term corresponds to the energy of electrons interacting with this self- 
organised magnetic field. The last term corresponds to the exchange energy and to fluctuations 
about the self-organised magnetic field (mean field) calculated to the order of e 2 . This term 
does not dominate over the two first terms in Eq. ( |119j ). Hence the fluctuations do not destroy 
the effects of the self-organised magnetic field. 



5 Ground state energy 

In the case T = the function H(r), which determines the self-organised magnetic field is 
known explicitly, as shown in subsection (3.4.). If the inequality (|68"D is in force, then the 
function H{r) is given by Eq. (|73|) . If the follwing inequality holds, 

H(0) > \£ F \ (121) 

then the function H(r) is given by Eq. (0). It should be noted that the latter case is ob- 
tained from the former in the formal limit k < — > 0, whence we continue our analysis using the 
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expression (|73|). Using this expression, or rather its Fourier transform ([73]), we can evaluate 
the energy E mag in Eq. ( |119| ) in a fairly explicit form. After straightforward calculations one 
obtains, 

E maa = At I —(kl + !> (122) 

3(4 + k%) - + k\) + 6(4 - k\f log ■ 



v \ 


r ^ 2 


2tt 2 1 


10m 


V \ 


fe 2 ■ 


2tt 2 1 


[8tt _ 



Before proceeding further, let us recall that the equations (|n]) and (^) have two positive 
roots k > and fc< if and only if the inequality (|80|) is in force. This inequality in turn implies the 
lower bound (^) for the Fermi energy Sp, which numerically is far above the Fermi energies 
for all known metals. 



The physically relevant case obtains when the inequality ( |121| ) holds. The corresponding 
formula for E mag can be obtained by passing to the formal limit k < —>■ in Eq. ( 122 ). Thus, if 
condition ( |121| ) holds true, then 

V \ h 2 3e 2 1 

E 0mag = hm E mag = — | _A£ - | (123) 

where the quantity k> is given either by the equation ([M]), which implies that A;> is in the 
interval (|92] ) or by one or the other of the equations (|99|) or ( |100| ), which restrict the quantity 
k~~, to be in one or the other of the intervals ( |102| ) or ( |103j ). 

Thus, for realistic metals, we expect the ground state energy to be given by Eq. ( |123| ), 
which originates from quasiparticles in the spin down states only, with energy spectrum £_(k). 

At this point it is appropriate to recall that we are supposed to calculate the energy per 
particle. The number of electrons is given as follows, 



N = E/d 3 x < ^l(x)Vax) > (124) 
d 3 k[/(£_(k)) + /(£ + (k))] 



V 

(2^)3 

where the functions f(£±) denote the appropriate Fermi functions ( |120| ). 



In the limit T — > the integral ( 124 ) can immediately be evaluated by use of the formula 
with the result, 

N = ^ k l + k <) (125) 

Thus the parameters k > and k < are related to the average particle density as indicated above. 
In particular, in the case k < — > we obtain the following relation from Eq. ( |125|) , 

^ = 6tt 2 4 3 (126) 
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Using the reult ( |126| ) in Eq. ( |123| ) we obtain the energy per particle, 

Eomag 3 2 f <^0 . 2 ^ 



e 2 _ fct, - — fc> (127) 
N 2 I 5 > 4tt > J v ; 

To obtain the genuine ground state energy, we still have to minimise the expression (127) with 
respect to the parameter fc>, which can vary within certain limits, as discussed above. One 
readily finds that the expression (|127|) has a minumum for the following value of k > , 



!5 1 1-875 ._ o . 

{k> )min = — = 128 

8 irao nay 



The corresponding minimum value of the energy per particle is, 



w \ / i \ 2 o 2 o 2 

^Omag \ _ K H « n oi/i e 



= -15 — — w -0.214 — (129) 

/mm V87T/ 2ao 2ao 



Numerically the result ( |129| ) corresponds to the following value, 

^Omag \ 



N A 



-2.91 eV" (130) 



The minimum ( |129| ) occurs at the value ( |128| ) for the parameter k > , which is within the 
limits in Eq. ( |103j ). Hence the minimum considered here corresponds to a negative value for 
the Fermi energy Sp, 

£ F <0 (131) 

In conclusion we wish to compare the results given above, with the ground state energy E 
of the unmagnetised case considered in Sec. 2. 

The ground state energy in question is obtained by adding the quantities E^ n and E exc h 
given by Eq.(|l4D and Eq. (|T5|), respectively, 

£ = £^/if£^ 5 i _±fM'lLf^(^l_^\ ( 132 ) 



2a I 5 V 4 J rl 2ir \ 4 / r s \ 2a n 1 r 2 

The parameter r s occuring in the equation ( |132j ) above, is related to the volume per electron, 
as detailed in Sec. 2, i.e. 

r, = - (133) 
a 

and 

4 Q V 

-7rr 3 n = — (134) 
3 iV V ; 

Thus, using Eq. Q126Q and the equations (|133|) and ( p.34| ) above, we obtain a relation between 
the parameters k > and r s , 

/97T\3 

r s fc>a = — (135) 
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Eliminating the parameter fc> from Eq. ( |127| ) by means of Eq. ( |135[ ) above, one obtains, 



e 2 AM3/97r\3l 3 /9vr\s 3 I e 2 N [ 3.508 1.731, 

e ^=^U{y) ^-^{t) ^r^br - ^ (136) 



The energy per particle of the magnetically ordered state is smaller than the corresponding 
energy per particle in the unmagnetised state, provided the parameter r s exceeds a certain limit 

<^^r s >r sc (137) 
N N 

A comparison of the two expressions ( |132|) and ( |136| ) above, yields a numerical estimate for the 
lower limit r sc in ( |137| ), 

r sc « 1.592 (138) 



The condition r s > 1.592 is satisfied for all metals, so that the energy per particle in the 
magnetically ordered state is smaller than the energy per particle in the unmagnetised state in 
all the relevant cases. 



Let us finally minimise the energy per particle (|132|) for the unmagnetised state with respect 
to r s . One obtains, 

N . 

The minimum value —2.91 eV for the energy per particle in the magnetically ordered state, 
which was obtained above in Eq (|129 ), is much closer to the experimental values for most metals 



-1.29 eV (139) 



[pjj| than the value —1.29 eV obtained above from the energy per particle in the unmagnetised 
case. The minimum in question occurs at 

8tt /9tt\s 

— — ^ 4.05 (140) 



as is readily found from Eqns. (128|) and (155) 



In fact, the binding energies of all metals exceed numerically considerably the value 1.29 eV, 
and are closer to the value 2.91 eV obtained for the magnetised case considered here. 



6 Summary and conclusions 

Starting form the simplest theory of metals, we have derived a particular mechanism by means 
of which the Coulomb interaction between electrons can generate a new kind of magnetically 
ordered state in a metal. The self-organised magnetic field geberated by this mechanism is 
represented by an order parameter which is a non-local field of a special kind. The self- 
organised field does not produce spin polarisations of the electrons, but generates particular 
long range correlations between electrons of opposite spin in different points x and x'. The 
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electron plasma in the model is described in terms of interacting quasiparticles of fermionic 
type. Each quasiparticle state is shown to be a special superposition of an electron spin up 
and spin down state, such that the probability of finding an electron with either spin up or 
spin down in each quasiparticle state equals |. Hence the spontaneously organised magnetic 
field does not generate polarisations of the electrons in the metal but only polarisations of the 
quasiparticles. 

The ground state of the electron plasma has been shown to be described by quasiparticles 
with a specific absolute spin-polarisation, but with no spin-polarisations of the electrons, for 
such values of e.g. the Fermi energy, which are valid for all existing metals. Hence the ground 
state in the simplest model of metals is not a ferromagnetic state as advocated by Bloch || , but 
rather a state involving the spontaneously generated non-local magnetic field derived in this 
paper. The ground state properties of this magnetised state give rise to theoretical numerical 
values for the average distance between electrons and for the binding energy which are in 
reasonable agreement with experimetal data. 

It ought to be possible to use our theoretical analysis as a starting point for experimental 
research in condensed matter physics aiming at observing and studying the magnetically ordered 
states predicted in this paper. 
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